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Abstract 



We investigate the relationship between ergodicity and asymptotic Gaussianity of isotropic 
spherical random fields, in the high-resolution (or high-frequency) limit. In particular, our 
results suggest that under a wide variety of circumstances the two conditions are equivalent, 
i.e. the sample angular power spectrum may converge to the population value if and only 
if the underlying field is asymptotically Gaussian, in the high frequency sense. These find- 
ings may shed some light on the role of Cosmic Variance in Cosmic Microwave Background 
(CMB) radiation data analysis. 
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^ ! 1 Introduction and background 

^ , 1.1 Overview 

The usual framework for proving asymptotic results in probability (for instance, central limit 
theorems or laws of large numbers) lies within the so-called Jarge sample paradigm, according 
to which more and more (independent or weakly dependent) random variables are generated, 
and the limiting behaviour of some functionals of these variables (e.g., averages or empirical 
moments) is studied. 

Physical applications, however, are prompting the development of a stochastic asymptotic 
theory of a rather different nature, where the indefinite repetition of a single experience is no 
longer available, and one relies instead on observations of the same (fixed) phenomenon with 
higher and higher degrees of resolution. 

One crucial instance of this situation appears when dealing with the statistical analysis of 
random fields indexed by compact manifolds, the quintessential example being provided by the 
case of the sphere S 2 . Indeed, we are especially concerned with issues arising from the analysis 
of the Cosmic Microwave Background (CMB) radiation, a theme which is currently at the core 
of physical and cosmological research, see for instance [12, 20J for textbook references, and 
[381 1391 121~] for further discussions around the latest experimental data. 



*We are grateful to Mirko D'Ovidio for some useful comments on an earlier version. 
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It is well-known that the CMB is a relic electromagnetic radiation providing a snapshot of the 
Universe at the so-called age of recombination, i.e. at the era when electrons in the primordial 
fluid arising from the Big Bang were captured by protons to form stable hydrogen atoms. Since 
the cross-section of hydrogen atoms is much smaller than for free electrons, after recombination 
photons can be viewed as diffusing freely across the Universe (to first order approximations). 
According to the latest experimental evidence, this has occurred some 3.7 x 10 5 years after the 
Big Bang, i.e. 13.7 billion years from the current epoch. Several experiments have been devoted 
to collecting extremely refined observations of the CMB, the leading role being played by the 



currently ongoing NASA mission WMAP (launched in 2001, see http://map.gsfc.nasa.gov/) 
and the ESA mission Planck, which is just now starting to operate after the launch on May 14, 
2009 (see |http : //www, sc iops . e sa. int/[ ). 

From a mathematical point of view, the CMB can be regarded as a single realization of 
an isotropic, zero-mean, finite variance spherical random field, for which the following spectral 
representation holds (see e.g. [lj or |22j) 

oo I 

T(x) = J2Y1 a imYi m {x) , x € S 2 . (1) 



1=0 m=-l 



Here, the collection 



{Y lm :l>0, m = -/,...,/} 



stands for the usual triangular array of spherical harmonics, which are well-known to provide 
a complete orthonormal system for the L 2 (S 2 ) space of square- integr able functions (with re- 
spect to Lebesgue measure) on the sphere - see [30] 03] 03] . In a loose sense, we can say 
that the frequency parameter I is related to a characteristic angular scale, say $i, according 
to the relationship ~ it /I. The (random) triangular array of spherical harmonic coefficients 
{aim '■ I > 0, m = —I, 1} are such that Eai m = and Eai m aii m i = Ci5\ 5™ , the bar denoting 
complex-conjugation and 5 h a indicating the Kronecker delta function. The non-negative sequence 
{Ci : I > 0} (not depending on m - see |29] as well as the forthcoming section) is the angular 
power spectrum of the spherical field (see for instance [21 [3]). 

As recalled above, our work deals with asymptotic issues, where the expression "asymptotic" 
has to be understood in the high-resolution (or high-frequency) sense. This means that we focus 
on the behaviour of the Fourier components 

l 

V (x) := a imYim (x), xe S 2 , I > 0, (2) 

m=—l 

associated with a fixed spherical field, as the frequency / grows larger and larger (plainly, each 
Ti is the projection of the field T into the orthogonal subspace of L 2 (5" 2 ) spanned by the 
spherical harmonics {Yi m : m = —I,..., I}). Note that this is the typical framework faced by 
experimentalists handling satellite missions as those mentioned above. Indeed, these missions 
are observing the same (unique) realization of our Universe on the so-called last scattering 
surface; more recent and more sophisticated experiments are then characterized by higher and 
higher frequencies (smaller and smaller scales) being observed. For instance, for the pioneering 
CMB mission COBE in 1989-1992 (which led to the Nobel Prize for Smoot and Mather in 2006) 
only frequencies in the order of a few dozens were recorded (i.e., scales of several degrees), a 
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limit which was raised to few hundreds by WMAP (i.e., approximately a quarter of degree) and 
is expected to grow to a few thousands with Planck (i.e., a few arcminutes). 

The principal goal of this paper is to enlighten some partial new connections between two 
high-resolution characterizations of spherical fields, that is, ergodicity and asymptotic Gaus- 
sianity. Roughly speaking (formal details are given in the forthcoming Sections 1.2 and 1.3), 
one says that the spherical field T is ergodic if the empirical version of the power spectrum of 
T (see formula (|3|) below) can be used as a consistent estimator of the sequence {C{\ (at least 
for high values of I). On the other hand, we say that T is asymptotically Gaussian, whenever 
suitably normalized versions of the frequency components of T\ exhibit Gaussian fluctuations 
for high values of I. As discussed below, these two notions are tightly connected whenever one 
deals with fields having an isotropic (or, equivalently, rotationally-invariant) law. 

Remark. For the rest of the paper, every random object is defined on a suitable (common) 
probability space P) . 

1.2 High-frequency ergodicity 

In what follows, we shall consider a real-valued random field T = {T (x) : x S S 2 } indexed 
by the sphere S 2 . The random field T satisfies the following basic assumptions: (i) the law 
of T is isotropic, that is, T has the same law as x ^ T (gx) for every rotation g £ SO (3) 
(here, we select the canonical action of SO (3) on 5 2 ); (ii) T is square-integrable and centered. 
Under assumptions (i)-(ii), the harmonic expansion ([1]) takes place, both in L 2 (P) (for fixed 
x) and in the product space L 2 (Q x S 2 ,P®d\), where A stands for the Lebesgue measure. 
Note that the last claim hinges on the fact that one can regard T as an application of the type 
T : Q x S 2 — > M : (uj, x) \— > T(u,x). As anticipated in the previous section, another useful 
property of T (easily deduced from isotropy - see e.g. [29]) is that the harmonic coefficients a; m 
are such that the power spectrum associated with T, defined as the collection {C; : I = 0,1, ...} 
(with C\ = E \ ai m \ 2 ), depends uniquely on the frequency indices I. 

In physical experiments (for instance, when measuring the CMB radiation), the power spec- 
trum of a given spherical field is usually unknown. For this reason, a key role is played by its 
empirical counterpart (called the empirical power spectrum - see for instance [H],[36]), which 
is given by 



An important issue to be addressed is therefore to establish conditions under which the distance 
between the quantities C\ and C\ converges to zero (in a sense that is defined below) when I — > oo, 
that is, when higher and higher frequencies of the expansion ([1]) are available to the observer. 
Although the asymptotic behaviour of spectrum estimators has been very deeply investigated for 
stochastic processes in Euclidean domains and under large sample asymptotics (see for instance 
[6l [22| [36] ) , only basic results are known in the high-resolution setting. 

For instance, it is immediate that the finite variance of T entails that, for every x £ S 2 , 




m=—l 



(3) 



ET {xf = 



(21 + 1) 



Q < oo, 



l>0 
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from which one deduces that Cj — * and also 

l>0 l>0 

By reasoning as in the proof of the Borel-Cantelli Lemma we therefore infer that, for any e > 0, 
P | lim supC/ > el < lim ^ P |q > e} < lim - J^C* = 0, (4) 

yielding in turn that both C\ and \C\ — Ci\ almost surely converge to zero as I —* oo. Plainly, 
since this result does not provide any information about the magnitude of the ratio \Ci — Ci\/Ci, 
it is virtually useless for statistical applications. In particular, one cannot conclude from (J3J) 
that the estimation of C\ based on C\ is consistent in a satisfactory statistical sense. 

Starting from these considerations, one sees that it is indeed necessary to focus on normalized 
quantities, such as the sequence 

m=— I 

Note that EC\ = 1, and also that the coefficient C\ is not observable (whereas Ci is). The 
sequence |q : I > ol can be used in order to meaningfully evaluate the asymptotic performance 

of any statistical procedure based on C\ . The following definition uses the coefficients C\ in order 
to define ergodicity. 

Definition 1 (HFE) Let T be an isotropic, finite variance spherical random field with angular 
power spectrum {Ci : / > 0} . We shall say that T is High- Frequency Ergodic (HFE - or 

ergodic in the high-frequency sense) if and only if 

lim E \di - l) 2 = lim E J § - 1 I =0. (6) 
Condition © implies of course that Ci = Ci/ Ci converges in probability towards the constant 

1. 

Remark. In some sense, the term "high-frequency consistency" seems to better describe prop- 
erty ([6]). However, in the statistical literature consistency is usually viewed as a property of a 
sequence of estimators, whereas here we deal with a property of the field T, so that we find 
the term ergodicity more suitable. Another way of formulating this point is to say that ([6]) 

characterizes the ergodicity of the "normalized empirical spectral measure" |C/ : I > o|, as I 
diverges. 

1.3 Ergodicity of Gaussian fields (and associated Gaussian fluctuations) 

As an illustration (and for future reference) we now test Definition [1] under the additional as- 
sumption that T is Gaussian. In this case, it is readily seen that, for every I > 1, the components 
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of the vector {ai m : m = are complex-valued and independent. Moreover, the random 

quantities aio/y/Cl, y/2Re (ai m ) / y/Ci and y/2lm (ai m ) / y/Cl (m = 0,...,l) are independent and 
identically distributed N(0, 1) random variables (these facts are well-known, see e.g. [2], [29] and 
the references therein). It is now easy to prove that 

m=— I 



in every norm LP, p > 1. Indeed, since af Q /Ci ™ x\ an d the set {2a 2 m /C; : m = 1,...,/} 
composed of i.i.d. x\ random variables independent of a/o (here, Xn denotes a standard chi- 
square distribution with n degrees of freedom), 



is 



r ~ i 2 l 

E\Ci-l\ = r^E 

I ' J (2/ + 1) 2 

2 



«2 ( ' \„ |2 

'0 _ 1 _!_ 9 J l Qhn l 

G 1 ^ C, 



21 + 1 l^oo 



0. 



and one can use the fact that, for polynomial functionals of a Gaussian field of fixed degree, all 
LP topologies coincide. 

We shall now provide (see the forthcoming Proposition [2]) a CLT that is naturally associated 
with the convergence described in ([7]). Note that, instead of using the classic Berry-Esseen results 
(see e.g. Feller [15]), we rather apply some recent estimates (proved in [32] and [33] by means 
of infinite-dimensional Gaussian analysis and the so-called "Stein's method" for probabilistic 
approximations) allowing to compare, for fixed I, the total variation distance between the law 
of the normalized random variable 



21 + 1 )Ci I /2/-M 



{Q-l}, 



and that of a standard Gaussian random variable. Recall that the total variation distance 
between the laws of two real- valued random variables X and Y is given by 

d TV (X, Y) = sup |P (X e A) - P (Y € A) | , 

A 

where the supremum runs over all Borel sets A. 

Proposition 2 Let N (0, 1) denote a centered standard Gaussian random variable. Then, for 
all I > we have 



^(i^ 1-1 , iV(0,l) <\It^tt> («) 



C, [ ' v ' ' ) - V 2Z + 1' 



so that, in particular, as I — > oo, 



21 + 1 



{a-l}^ N(0,1). (9) 
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Proof. We have 




+ ^ 2 ( ReQ ^} 2 + {Ima; m } 2 _ ^ + ^ 



where {xi m } are a triangular array of i.i.d. standard Gaussian random variables. Standard 
calculations yield that 



CU7U4 



V2 



21+1 ( T i _n 



m=l 



12 



21 + 1' 



where cumj stands for the jth cumulant. Now recall that in [33] it is proved that, for every zero 
mean and unit variance random variable F q that belongs to the qth Wiener chaos associated 
with some Gaussian field (q > 2), the following inequality holds: 



( / /r (7 v .Yi(). h) < 2\ . 



The result now follows immediately, since each variable ^/ { %[ ^ as un ^ var i ance > an d 
is precisely an element of the second Wiener chaos associated with T. ■ 

It is simple to verify numerically that the convergence ([9]) takes place rather fast. For 
instance, for / = 100, the bound in total variation is of the order of 2%, while for / = 1000 we 
deduce an order of 0.6%. 

We stress that the previous results heavily rely on the Gaussian assumption, and cannot 
be easily extended to the framework of non-Gaussian and isotropic spherical fields. The main 
reason supporting this claim is contained in the references [21 [3], where it is shown that, under 
isotropy, the coefficients a\ m are independent if and only if the underlying field is Gaussian, 
and this despite the fact that they are always uncorrelated by construction. In other words, 
sampling independent, non-Gaussian random coefficients to generate maps according to ([T|) will 
always yield an anisotropic random field. The dependence structure among the coefficients 
{aim} is i n general quite complicated, albeit it can be neatly characterized in terms of the group 
representation properties of SO{2>) (see [28] and [29]). In view of this, to derive any asymptotic 
result for C\ under non-Gaussianity assumptions for T, is by no means trivial and still almost 
completely open for research. 

1.4 High-frequency Gaussianity 

A different form of asymptotic theory has been addressed in an apparently unrelated stream of 
research, for instance in |28j. 

Definition 3 (HFG) Let T{x) be an isotropic, finite variance spherical random field, and recall 
the notation (CP and (dp. We say that T(x) is high-frequency Gaussian (HFG) whenever 

Tl(x) law 



^Var{Ti{x)} 

for every fixed x £ S 2 . 



N(0, 1) , as I oo, (10) 
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Remark. It is more delicate to define HFG involving convergence in the sense of finite di- 
mensional distributions. Indeed, in [28] it is shown that, even if relation (jlOp holds, the finite- 
dimensional distributions of order > 2 of the field x *— > T[(x)/ \/Var may not converge 
to any limit. 

It is clear that a Gaussian field is asymptotically Gaussian: however, as shown in [28J, char- 
acterizing non-Gaussian fields that are HFG can be a difficult task, even if the underlying field 
T is a simple transformation (for instance, the square) of some Gaussian random function. Con- 
ditions for the HFG property to hold in some non-Gaussian circumstances are given in [28], by 
using group representations - yielding some interesting connection with random walks on hyper- 
groups associated with the power spectrum of T. We stress that the possible existence of HFG 
behaviour entails deep consequences on CMB data analysis. On one hand, in fact, parameter 
estimation on CMB data is largely dominated by likelihood approaches, whence an asymptot- 
ically Gaussian behaviour would great simplify the implementation of optimal procedures. On 
the other hand, testing for non-Gaussianity is a key ingredient in the validation of the so-called 
inflationary scenarios, and the possible existence of high frequency Gaussian components for 
non-Gaussian models might set a theoretical limit to the investigation in this area. 

1.5 Purpose and plan 

Our purpose in this paper is to investigate the relationships between the HFG and HFE prop- 
erties under an assumption of Gaussian subordination, that is, by considering fields T that can 
be written as a deterministic function of some isotropic, real-valued Gaussian field. We will 
mainly focus on the case of polynomial subordinations, where the polynomials are of the Her- 
mite type. Note also that Gaussian subordination is the favoured framework for CMB modeling 
in a non-Gaussian setting (see e.g. [4], [17], [45]). 

Our main finding is that, despite their apparent independence, the HFG and HFE properties 
will turn out to be very close in a broad class of circumstances, suggesting that ergodicity (and 
hence the possibility to draw asymptotically justifiable statistical inferences) and asymptotic 
Gaussianity are very tightly related in a high-resolution setting. This may lead, we believe, to 
important characterizations of Gaussian random fields, and to a better understanding of the 
conditions for the validity of statistical inference procedures based on observations drawn from 
a unique realizations of a compactly supported random field, as in the spherical case. 

The plan of this paper is as follows: in Section 2 we state and prove our main result, 
establishing necessary and sufficient conditions for ergodicity and Gaussianity and exploring the 
link between them. Indeed, these conditions turn out to be extremely close, so that in Section 3 
we can indeed discuss more thoroughly a special case of practical relevance, namely the quadratic 
case. Section 4 is devoted to further discussion and directions for further research. 

2 A general statement about Gaussian subordinated fields 

The two notations ([I]) and (|2]) are adopted throughout the sequel. Let us first recall a few basic 
facts and definitions. 

(I) The first point concerns a characterization of isotropy in terms of angular power spectra. 
Indeed, as discussed in [13129], if a random field is isotropic with finite fourth-order moment, 
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then there exists necessarily an array |7^ 4 (L)| such that 

cum {cii irni , CLi 2Tn2 , Q>l 3ms , 0/4^714} (11) 

= ^ (_1)M ( L\ m 2 M ) ( m 3 m 4 — M ) ( 2L + ^fa'W • 



In general, the symbol cum {Xi, X m } denotes the joint cumulant of the random variables 
X\, X m . Also, we label as usual |7^ 4 (L)| the cumulant trispectrum of the random field 

(see for instance [IT], [29]); as made clear by our notation, the quantity 7^' 3 / 4 (L) does not depend 
on mi, m-2, m.3, 7714 (this phenomenon is analogous to the fact that the power spectrum only 
depends on the frequency I - see |29] for a discussion of this point). On the right-hand side 
of (jlip . we have also introduced the well-known Wigner's coefficients, which arise in the group 
representation theory of SO (3) and are discussed at length in many excellent monographs on 
the quantum theory of angular momentum - see e.g. [23, 43, ;43]. For future reference, we also 
recall that the Wigner's 3j coefficients are equivalent, up to normalization and phase factor, to 
the Clebsch-Gordan coefficients given by (see e.g. 



= (-1)"—' V5E + l(£ £ L M ) • (12) 

As noted by [17] . geometrically the multipoles (h,l2,hih) can De viewed as the sides of a 
quadrilateral, and L as one of its main diagonals; L is also the shared size of the two triangle 
formed by the corresponding pairs of side. Clebsch-Gordan coefficients ensure that the triangle 
conditions are satisfied, indeed they are different from zero only if. l\ < I2 + L, I2 < l\ + L, and 
L < h + l 2 . 

(II) We shall sometimes label a point x of the sphere S 2 in terms of its spherical coordinates, 
that is, x = 99), where < 1? < ir and < (p < 2ir. 

(III) Easy considerations yield the important fact that, for any isotropic random field T, 



Ti(ti,tp) la =? Ti(N) =J2"l™Yi m (N) ^ a l0] j ^±1 , 

Ira 

where we denote by N := (0, 0) the North Pole of the sphere and by " l == " the equality in law 
between two random elements. 

(IV) It is immediate that, if T is isotropic, then for every deterministic function F the sub- 
ordinated random application x 1— > F (T (x)) is also isotropic. Moreover, if F(T (x)) is square 
integrable, then F(T(-)) also admits a harmonic expansion analogous to ([T]). One specific in- 
stance of this situation is obtained by choosing T to be Gaussian and isotropic, and F to be any 
of the Hermite polynomials {H q : q > 0} (in this case, one talks about a Gaussian subordination 
of the Hermite type). We recall that the polynomials H q are such that H q = 5 q l, where 1 stands 
for the function which is constantly equal to one, 8° is the identity, and 5 q (q > 1) represents the 
qth. iteration of the divergence operator 5, acting on smooth functions as 5f (x) = xf (x) — f (x) . 
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For instance Hq = 1, H\ (x) = x, H2 (x) = x 2 — 1, and so on. When T is Gaussian, we adopt 
the notation 

00 

H q {T(x)):=T q (x) = Y J T hq (x), x G S\ q > 2 , (13) 



where 



1=0 
1 



T l;q( X ) = ^2 a lm-qYl m {x) (14) 
m=—l 

is the Zth frequency component of T q , with ai m - q the associated harmonic coefficients. We 
shall also write {Cz ;g : I > 0} and {^(L;g)}, respectively, for the power spectrum and for the 
cumulant trispectrum of T q , as introduced at Point (I). According to [281 Theorem 3], one has 
that Ci- q admits the following expansion in terms of the power spectrum {C{\ of T: 

cw t o.-AaSrfn^} E R£w^f- (15) 

h,...,J g =0 U=l J Li...L 9 _ 2 

where C^ 1 ^ 2 '.^'^^ 2 ' 1,0 indicates a convolution of Clebsch-Gordan coefficients, that is 

Ai Ap_2 

^-,Ai,A2,...,A p _i;/j V ■> \ ^-,Ai,/ii fj\'2,H2 _ _ _ /--(Ap_i,/i flfi) 

fi\=— Ai Mp-2 =— A p _2 

(see |28| and [29J for more details on these convolutions, which can also be viewed as probability 
amplitudes in alternative coupling schemes for quantum angular momenta, compare [5]). 

(V) An easy but important remark is the following. Since the expansion ([1]) is in order, the law of 
a centered isotropic Gaussian field T is completely encoded by the power spectrum {C; : I > 0}. 
This is a consequence of the fact that, in this case, the array {a; m : I > 0, m = 0,...,l} is 
composed of independent Gaussian random variables such that: (i) a/o is real-valued, and (ii) 
for every m > 1, the coefficient a; m has independent and equidistributed real and imaginary 
parts. 

As anticipated, we shall now prove some new connections between HFE and HFG spherical 
fields (see Definitions [1] and [3j) , in the special case of fields of the type T q , as defined in (I13p . 
In particular, our main finding (as stated in Theorem H]) Note that the conditions appearing in 
the following statement involve the coefficients Ci- q given in (|15p . and that these coefficients are 
completely determined by the power spectrum of the underlying Gaussian field T. 



Theorem 4 Let q > 2, and define T q according to ( G2P , where T is Gaussian and isotropic. Let 
ylsh 



T^f 4 (L; q) be the reduced trispectrum of T . Introduce the notation 



w u (L) := (Ct %) 2 and w 2l (L) = 

in such a way that 

21 21 
L=0 L=0 

Then, the following holds. 
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1. the random field T q is high-frequency Gaussian if and only if 



limXSW^£^ = 0. (17) 

l=o 

On i/ie oi/ier hand, T q is high-frequency ergodic if and only if 



Hmf>^)^^ = 0. (18) 



I— »oo ^— ' C; 

l=o 



Before proving Theorem HI we shall note that {C^ } 2 is different from zero only for L even, 
and T,f(L) is not in general positive- valued. Moreover, in view of the forthcoming Lemma 



also in (|18|) the sum runs only over even values of L. 
Lemma 5 T}}(L) is zero when L is odd 

Proof. From Eq. (17)], we infer that, in general, 

T^(L) = (-l) l ^ +L T^(L) . 
Considering the case h = h = h = h = h we obtain the desired result. ■ 

Proof of Theorem [4] (Proof of 1.) Consider the random spherical field 

(<?,¥>) ~T, ; ,(0,p) := SlgS^L :, (#,</>) G [0,vr] x [0,2vr), 
^ar{T, ; ,(JV)} 

where iV is the North Pole, and observe that, by isotropy and for every (#,<£>), 

rf, / o \ law OlQ 



The field Ti ;q is mean-zero and has unit variance: since it also belong to the qth Wiener chaos 
associated with T, we can deduce from the results in [31] it is asymptotically Gaussian if and 
only if 

lim -p^-cum A {a m . q } = . 



l^oo C lq 



As discussed e.g. in |17| and [29], isotropy entails that we can write the fourth-order cumulant 

as 

cWW = B-Wo o m)(o o _ L M )w< + Wti(L) 
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so that the field is asymptotically Gaussian if and only if 

2 



Since relation (|12p is in order, we write 

I I L x 




(2L + 1) = {C$ } 



entailing in turn that 

2 2/ 2/ 



E(o o o) ^ + 1 ) = EW 2 = E E («} 2 -i 

r V 7 L=0 L=0M=-L 



where the second equality follows from the fact that Clebsch-Gordan coefficients Clj^^ are 
different from zero only for 771,3 = m i + m 2, and the third equality is a consequence from the 
orthonormality properties of the coefficients (which are the elements of unitary matrices whose 
rows are indexed by 7771,7772 and whose columns are indexed by £3,7713). We therefore have 



^-cum 4 {a m . g } = -4- E {CWW 



yielding the desired conclusion. 

(Proof of 2.) On the other hand, we obtain also 



1 1 r - - 2 1 \^Jp| l 2 \ 2 

(2/ + l) 2 C^~~ cum { aimi '<Q ,aimi ''y> aim2 '' ( i ,airn2 '' c if ^ (21 I p ? ^ 1 l a ^;gl J v^U) 

^ ' raim2 ^ ' J)9 m 

^ ^ 2 

= (2/ + l) 2 C 2 "" ^ ' ( _ ^) 1+ 2 CUm {a; mi;g , CLl ) -mi;q, a lrri2;qi a l,-rri2;q} + 7^ ~ ^\ (21) 
^ ' mirri2 ** 



■ii ■ \r 2 cr 2 ^ 11 v ///, »/■» a/ 

Z Z L \ . ^ W/ ^ . 2 



x 



x »r ](2L + 1)T, ,{L) + — (22) 

» -7771 -7T7 2 ~M J ^ ' 11 V ' (21 + 1) V ^ 

2Z 

piw4x?o (2L + 1)7I?<i)+ (^)' <23> 

L even 



It is simple to notice that 

2/ 



<T(2L + 1) = 2 2l{2l + 1) + 21 + 1 = (21 + 1) 



L=0 
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so that we have 



E 



_ 1 I = 2 £ w l l L Tll(L) + — ^— , where m > and £ m = 1. 

l ' q J L=0 ^ ' L=0 

L even L even 



The result now follows immediately. 
Remark. Note that 



L0 2 _ (2L + 1)(^!) 2 (L!) 2 (2?-L)! 1 



|!)2 (2/ + L + 1)! " (2L + 1) 



(2L + 1) 1 



(2Z + l) 2 - 21 + 1 

Note also that in the Gaussian case (e.g., q = 1) we have T^(L) = 0, whence 

2 

2 

, 




{21 + 1) 



as expected. 



The previous result strongly suggests that the conditions for asymptotic Gaussianity (HFG) 
and for ergodicity (HFE) should be tightly related. Indeed we conjecture that HFE and HFG 
are equivalent in the case of Hermite type Gaussian subordinations (and most probably even 
in more general circumstances). However, proving this claim seems analytically too demanding 
at this stage, so that for the rest of the paper we content ourselves with a detailed analysis of 
quadratic Gaussian subordinations. In particular, we believe that the content of the forthcoming 
Section [12] (which is already quite technical) may provide the seed for a complete understanding 
of the HFG-HFE connection. 



Remark 6 It should be noted that the reduced trispectrum satisfies (see Eq. (16)]) 

7if(L') = B 2L + 1 ){ \ \ V 

In the previous remark, we introduced the well-known Wigner's 6j coefficients, which inter- 
twine alternative coupling schemes of three quantum angular momenta (see [5], |43j for futher 
properties and much more discussion). Their relationship with Wigner's 3j coefficients is pro- 
vided by the identity 

(a b e 1 = ^ ,_ Y+f+e+(f) ( a b e \ / c d e \ ( a d f \ ( c b f \ 
led//- ^ K > \a (3 e )\ 1 8 -e J\a 5 -<P ) \ <y [5 </> J 

(24) 

(see [33], Chapter 9, for analytic expressions and a full set of properties). 
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3 The quadratic case 

3.1 The class 3D and main results 

As anticipated, the purpose of this section is to provide a more detailed and explicit analysis of 
the quadratic case q = 2. For simplicity, in the sequel we consider a centered Gaussian isotropic 
spherical field T such that Var(T(x)) = + l)Cj/47r = 1, where {Cz} is as before the power 

spectrum of T. We start by recalling the notation 

oo 

T 2 {x) = H 2 (T(x)) = ^2 ^2 a hmi a hm2 Y hmi (x)Yi 2m2 (x) - 1 , (25) 

h,l 2 =l rai»2 

where T is isotropic, centered and Gaussian. Our first result can be seen as a consequence of 
formula (|15p (or, more generally, of the results of |28j). Here, we provide a proof for the sake of 
completeness. 



Lemma 7 The angular power spectrum of the squared random field $25\i is given by 



C/ ; 2 — E\ai m;2 \ 2 — 2 C; t C, 



h 



hh 



h l 2 I \ (2/ + l)(2/ + l) 







47T 



Proof. Recall first that Yqq(x) = (4ir) 1//2 , see |43| . equation 5.13.1.1. Hence, in view of 
25]), we have that, for / = 0, 



^00:2 



s- 



^ ^ a hm 1 ai 2 m2 Y hm 1 (x)Yi 2m2 (x) - 1 ^ Y 00 (x)da 

lxl 2 m\m2 



\=yZ ^2 a hm 1 ai2m 2 l I Y hmi {x)Yi 2rrl2 (x)dx\ - y/A-K 

k=J2^2 a hra 1 ai 2m2 5\ 2 1 5 r ^ 2 1 ~ 
hh 

-y 



l x l 2 m\m 2 



47T . 



Im 



It follows that 



and 



Im K Im J 



EH 2 (T(x)) = Ej2a lm . 2 Y lm (x) = Ea 0{) , 2 Y m (x) = 



z=o 



the second step following because Eai m = for all / > under isotropy (see [2]). Indeed we 
have (from (|25p . and in view of (|14p ) 



;2 = / ^2 a iimi a «2m 2 ^ mi (x)yz m2 (2;)y im (x)(ia 



i 1 Z 2 mim2 
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v- I h h I \f h k l\ / (2/ 1 + l)(2/ 2 + l)(2/ + l) 
A. A. a ^ a ^ ( mi m2 -m J i J V 4^ " (26) 



l lt l 2 =l m 1 m 2 

Note that the constant term —1 has no effect for I > 1, because 



/ Y lm (x)dx = for all I > 1 . 
Js 2 



Now 



h h I \(hh l \ /(2/i + l)(2/ 2 + l)(2/ + l) 



/ V 4vr 



h ^ mi ^ 



/i Zi £ 
m-i — mi 



ii 

in view of the well-known properties ([43, Eq. (8.5.1.1) and Eq. (8.7.1.2)]) 



h h 




) = , 1 , V(-i) mi ( h h 11 = V 7 ^!^ • 



Hence we have, as expected, Eai m -2 = for all I. Furthermore 

h h I \ fh h l\ /(2/i + l)(2/ 2 + l)(2/ + l) 



E\a lm , 2 \ 2 = eIJ2H ^^(^ ^ J m )(o o o) 



Z1Z2 Wl 7712 





■ h 


)( 





I'l 


'2 










to 




; ■ 



47T 



- - (A 1' 2 I \ { l[ % 1 \ / (2? / 1 + l)(2/ 2 + l)(2/ + l) 
2^ 2^ a 'imi a '>2 I m > m > - m ) I ) V 4tt 

l[l' 2 m' x m' 2 ' ' 

I mi m2 - m J I 7 4vr 

/1/2 mim2 v ' v ' 

r f Zl /2 z A 2 (2^i + l)(2/ 2 + l) 

and the proof is completed. ■ 
Remark. Note that 

Var{T\ X )} = Y J 2 -^C l 

2/ + 1 / h l 2 



4?r 



3s w 2ii+ r 2+i) fe 

ilia " I I 

nV^^ « (2/l + 1)(2Z 2 + 1) r T/ [rr , U1 2 

2 2^Q 1 C/ 2 = 2[yar{T(x)}J , 
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as expected from standard property of Gaussian variables. Here we have used again 

B 2 < + i)(o o J)' = 1 - 

Our strategy is now the following. We shall first define a very general class, noted D, of 
quadratic models in terms of the power spectrum of the underlying Gaussian field, and then we 
shall show that the two notions of HFG and HFE coincide within 2) . 

Definition 8 The centered Gaussian isotropic field T is said to belong to the class D if there 
exist real numbers a, (3 such that 

1. a£R and f3 > 

& E£o rQ+le_/3/ < 00 

3. there exists constants c\,C2 > such that 

C Ci 
< ci < lim inf m < lim sup < C2 < oo (27) 

Remarks. (1) As a first approximation, the class D contains virtually all models that are 
relevant for CMB modeling in the case of a quadratic Gaussian subordination. For instance, 
Sachs- Wolfe models with the so-called Bardeen's potential entail a polynomial decay of the C[ 
(/3 = 0), whereas the so-called Silk damping effect entails an exponential decay of the power 
spectrum of primary CMB anisotropics at higher /. We refer again to textbooks such as |12j . 
[13] for more discussion on these points. 

(2) Note that Condition 2 in the definition of £> implies that the parameters a, [3 must be 
such that either (3 = and a > 2, or (3 > and a £ M (with no restrictions). 

The next statement is the main achievement of this section. It shows in particular, that the 
HFG and HFE exhibit the same phase transition within the class D 

Theorem 9 Let T2 = H2 (T), where the centered Gaussian isotropic field T is an element of 
the class 2). Then, the following three conditions are equivalent 

(i) T 2 is HFG 

(ii) T 2 is HFE 

(iii) (3 > and a £ R. 
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3.2 Proof of Theorem M 

From [281 Section 6] , we already know that Conditions (i) and (iii) in the statement of Theorem 
[9] are equivalent. The proof of the remaining implication (ii) (iii) is divided in several steps. 

We start by showing that, if (iii) is not verified, then the angular power spectrum of the 
transformed field, under broad conditions, exhibits the same behaviour as the angular power 
spectrum of the subordinating field. 

Lemma 10 Suppose (3 = and a > 2, then 

^-cA < d, 2 < -i- {2((a - 1) + C(a)} C l/2 = 0(Q) , 

4-7T C 2 C\TT 

where £(.) denotes the Riemann zeta function. 
Proof. We have 

2 



r (h h I X (2h + l)(2l 2 + l) ^ f h i 2 l V (2Zi + l)(2l a + l) 

Z^ C ^ I J 4^ " 2 ^ ° lA I J 4^ 

hh 7 h<h V 7 

f h h I V (2*1 + l)(2/ 2 + 1) 

< 2-C l/2 ^ CU 1 — 

1 /i<i 2 v 7 



4vr 

because (Zi V l 2 ) > 1/2 by the triangle conditions and sup; 2 >// 2 Ci 2 /Ci/ 2 < c 2 /c\. Now 

h l 2 I V (2Zi + l)(2Z 2 + l) 



E ( '„* '„ 2 i ) 



47T 



(2Z 1 + 1) ^ ( h l 2 iy ^ ( 2 f 1 + l) 



h h 

More precisely 



47T 



£ c, ^ £ £ E(» + 1 ' r ° £ £ < 2 C(» - D + «a)} 

I 



Hence 



C, ;2 <-^-{2C(a-l)+C(a)}Ci /a . 



The upper bound is then established. For the lower bound, it is sufficient to show that 

2 " f ' h l 2 I \ 2 3(2Z 2 + 1) 



/ 4vr 



I o o o J ^ >2^ c i^ 

ilia V 7 h 

> o 9 a r cf^/Z 1 Z 2 Z \ 2 (2Z 2 + 1) 3x2° 



h 

as claimed. ■ 
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Loosely, the previous Lemma [TOl states that, under algebraic decay, the rate of convergence to 
zero of the angular power spectrum is not affected by a quadratic transformation, i.e. Ci-2 — C\. 
The following result holds for fixed I, and it is therefore not related to the high-frequency 
asymptotic behaviour of the power spectrum {C{\ (see [28J for related computations). Note 
that we use the notation 

I 



2 n C W 

«Zm;2| 

m=— I 



Lemma 11 Let T2 be defined by \25)) . Then we have 

£{G ;2 -i} 2 



16 X^r*n r f 1 ^ ^ I \V h h I ^ 2 (2Zx + 1)(2Z 2 + 1)(2Z 3 + 1) . 
Cf 2 ^ L h L h^h [ 00 )[ 00 ) (4^2 + K W > 



i;2 i x i 2 i z 
where for all I = 1,2, .... 



< R(l) < ' 



21 + 1 



Proof. In the sequel, we shall use repeatedly the unitary properties of Clebsch-Gordan 
coefficients, i.e. 

V ( 1 1 L \( 1 1 L ' \ = 6 1^_ (28) 
^ I mi m 2 M I \ mi m 2 M' ) 2L + 1 

Recalling [201 Ell we need to evaluate 



(21 + 1) 2 C 2 ^ > cum {a/ mi , ai mi , a,i m2 ai m2 } 



(21 + iyq 2 



^2 (-l) mi+m2 cum {a lmi , ai _ mi , a /m2 oj _ m J 



mim2 



Now 



CUm {cLlrrii 5 1; _ mi , a,i m2 ai _ m2 } 



cum < 



W2 /ilM2 



Zi Z 2 I Uli lj H / (2Zi + l)(2f 2 + l)(2Z + l) 
/!| "^ 2 l/n /u mi OJV 1- 



v- v- /l 3 ! 4 I \ ( h k I \ / (2Z 3 + l)(2/ 4 + l)(2/ + l) 

2^ Z> a '3w^4M4 1 M o J y 47r 

1 1 on a \ ' ' / \ / 



^2 ai ^o a ^6 
W8 



h k I \ ( h k l\ l(2k + l)(2k + l)(2l + l 



H5 fiQ 7TT-2 J V 1 V 47T 



h h I \( h h l\ / (2/ 7 + l)(2/ 8 + l)(2/ + l) 
H7 us ~m 2 J I J V 4tt 
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and counting equivalent permutations 

, , , x / , , , x (2/ + l) 2 n (2li + 1) 

'1 h ' W ll (2 I \ i=l 



E E ( _ 1)w+ „ + „ + „ CIIC1JclJcll ( fa fa ^ ) ( fa fa I ) 



(47T) 



*1 ^3 * W il ^3 M / h I \f k h l\ ( k h I \( k h 

-fil -(13 -mi J \ J I (14 (13 m 2 J \ J \ -(14 -(12 m 2 J \ 



+8 E E (-ir + « + <" + <"«,c lA q, ( <> £ „' ) ( 5 J J ) 



+ 1) 2 n (2i< + 1) 

i=l 



(47T) 



(13 (*4 -mi / 1 J \ -(ii -(i 3 m 2 y 1 J \ -(14 -(i 2 -m 2 j \ 



+ 8 E E (-ir i+ ^ a+ ^c h c l2 c h c u ( h h 1 ) 

h k l\( h h I \( h h I 
) \ -m (13 -mi j\0 

4 

( 2i + !) 2 n + 1) 

Z2 44 I \ I k k I \ 1 h k I \ ( k k I \ i=i 



-(12 (14 m 2 J \ J \ —(13 —(14 —m 2 J \ J (4tt) 2 
=: 8{A(mi, -mi,m 2 , -m 2 ) + B(mi, -mi,m 2 , -m 2 ) + C(mi, -mi,m 2 , -m 2 )} . 
For the first term, note first that (— l) m i+ m 2+Mi+M2+M3+M4 = ^ because the exponent is neces- 



sarily even by the properties of Wigner's coefficients. Moreover, applying iteratively ([28 

£ A(mi,-mi,rn 2 ,-m a ) 



mitri2 



E E E c h c h c h ci 4 ( q z 2 j 



, (2/ • I! 2 lli2/, • I) 

tl <2 t \ / £l 13 * \ i=l 



y (4vr) 5 
Z4 Z3 I \ ( k I3 I \ ( k k I \ f k k I \ ^3^il 



hhhU m 2t*2 M3M4 

a*4 ^3 m 2 y ^ y v, — ^4 — ^2 m 2 y ^ o o o j 2/3 + 1 

n n^n ( l ^ ^ iVf k k I \ 2 (2Zx + 1)(2Z 2 + 1)(2Z 4 + 1)(2Z + l) 2 

- 2^ c h^u ( I I o o J 7^ • 

Likewise, for the second term we note that (— 1)^1+^2+^3+^4 = \^ anc [ us i n g (|24p 

(-l) mi+m2 B(mi,-mi,m2,-m 2 ) 

m\m,2 
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^2 C h Ci 2 C h C h 

I1I2I3I4 



h h I 1 ( h h I 



Z3 U I 




h h 1 
000 



1 4 I2 1 



I4 I2 1 
000 



000 



{21 + if n (2«< + 1) 

i=l 



(4vr) 2 



Now by Cauchy-Schwartz inequality and recalling that 



h h I 
l 2 I4 I 



< 



1 



21 + 1 



for all hjfa, h,h 



the previous quantity can be bounded by 



21 



< 



1 r> n n r> ( ^ ' 2 ^ \ ( ^ ^ ^ \ f ^ ^ 3 ' 1 

— 2^ c^ 2 o^ 4 1 Q J I 0M0 J 



c h c h c h c h 

hfateU 

^2 C h Ci 2 C h C h 



I4 I2 I 




(21 + If n (2ii + 1) 
i=l 



(47T) 



il / 2 I 





h h I 





l 3 U I V (2Z 3 + 1)(2Z 4 + 1) (2h + l)(2/ 2 + 1)(2Z + 1) 




47T 



4tt 



l 2 l 4 I \ ( 2/ 3 + 1 )( 2 ' 4 + X ) ( 2Zl + l ^ 2h + 1 )( 2Z + !) 




47T 



4tt 



1/2 



1/2 



2^ + 1^2 



whence 



{2i + \yc? 2 S 5(mi ' 



-mi,m 2 , — m 2 ) 



< 



2/ + 1 



It is easy to see that ^ mim2 -mi,m 2 , -m 2 ) = E mim2 C(mi, -m 1 ,m 2 , -m 2 ). In view of 

[20l [21~1 the statement of the lemma follows easily. ■ 

The proof of Theorem [9] is now concluded by the following lemma. 

Lemma 12 If (3 = and a > 2, i/ien 



lim inf E \ C l;2 - 1 > > C 

l— >oo 



{Ci;2 " l} 
J/ (3 > and a is real, £/ien 



" 2 {2C(a-l)+C(a)}2 a J> > 0. 



2cfvr 



Hit 1 E{c Z;2 -l} =0 . 



Proof. For the first part, from Lemma [TT] we can focus on 



1 n 2 n n ( h h I \ 2 ( h h I \ 
Cpl ^ G h G b c ia I M J 



h k l \ (2/i + l)(2/ 2 + l)(2/ 3 + l) 
/ (4tt) 2 
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' £(2 il + i)(2 i2 + i,c llC , 2 (<> o o) 2 S Q A(o o o) W 



«;2 w 2 

which is larger than 



da 



2/ 3 i V (2/3 + 1) 



/ (4tt) 2 

-~ch~V 2 } lo oj Mo oj -ii^r = -^r 



'J;2 h \"""/, s \""-/ V-'V ^Z ;2 

Now we have proved earlier that in the polynomial case, C/ ; 2 — Ci ~ Z~ a , so the previous ratio 
does not converge to zero and Cj. 2 cannot be ergodic; the lower bound provided in the statement 
of the Lemma follows from previous computations and easy manipulations. 
For the second part of the statement, it is sufficient to note that 

2 ' h h l\ 2 (2/3 + I) 



^L^(2/ 1 + l)(2/ 2 + l)C /l Q 2 ( h Q h Q J) 52C h Q 



C 2 A^v— 1 1 / V / (4vr) 2 

< sup^^ + 1)- 1 Ez 3 r tl r t , {cj% h0 } 2 ^ BUPtl e,, r tl r t3 KVq} 2 ^ 

SilJs r 'i r ^3 {C/iO/30} Shis r 'i r '3 {^liOJ 3 o} 

so the condition is met, just as for the standard case of [28]. ■ 

Remarks. (1) By inspection of the previous proof, we note that we have shown how the 
sufficient condition for asymptotic Gaussianity (HFG) is also such for ergodicity (HFE). More 
precisely, we have proved that 

hm sup „ = hm sup-P(Zi = h\Z2 = fa) = , 

where {Z{\ is the Markov chain defined in [28, Eq. (57) and (58)] is a sufficient condition for the 
HFG (see [281 Proposition 9]) and also a sufficient condition to have lirn^oo E iC\ — 1 j =0. 

(2) In principle, the case q = 3 can be dealt along similar lines. 

(3) (On Cosmic Variance) Loosely speaking, the epistemological status of Cosmological 
research has always been the object of some debate, as in some sense we are dealing with a 
science based on a single observation (our observed Universe). In the CMB community, this 
issue has been somewhat rephrased in terms of so-called Cosmic Variance - i.e., it taken as 
common knowledge that parameters relating only to lower multipoles (such as the value of 
Ci, for small values of I) are inevitably affected by an intrinsic uncertainty which cannot be 
eliminated (the variability due to the peculiar realization of the random field that we are able to 
observe), whereas this effect is taken to disappear at higher I (implicitly assuming that something 
like the HFE should always hold). Our result seem to point out, apparently for the first time, 
the very profound role that the assumption of Gaussianity may play in this environment. In 
particular, for general non-Gaussian fields there is no guarantee that angular power spectra and 
related parameters can be consistently estimated, even at high multipoles - i.e., the Cosmic 
Variance does not decrease at high frequencies for general non-Gaussian models. 
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4 Discussion and directions for further research 



This paper leaves many directions open for further research. We believe the results of the 
previous two sections point out a very strong connection between conditions for High Frequency 
Ergodicity (HFE) and High Frequency Gaussianity (HFG) for isotropic spherical random fields. 
It is natural to suggest that equivalence may hold for Gaussian subordinated fields of any order 
q, or even more broadly for general Gaussian subordinated fields on homogeneous spaces of 
compact groups. Indeed, in this broader framework it is shown in [3] that independence of 
Fourier coefficients implies Gaussianity, which is the heuristic rationale behind our results here. 

The connection between the HFE property can also be studied under a different environment 
than Gaussian subordination. Consider for instance the class of completely random spherical 
fields, which was recently introduced in [8j [9]. Following the definition therein, we shall say 
that a spherical random field is completely random if for each I we have that the vector a/. = 
(a^_;, ...,au) is invariant with respect to the action of all matrices belonging to SU{21 + 1) and 
verifies a\ m = (— l) m a/ m . Because of this, the vector a/, is clearly uniformly distributed on the 
manifold of random diameter ^ \aim\ 2 = Ci, or equivalently, introducing the (21 + 1) vector Ui 

^ _ 1 I \/2Re ag \Z2Re a t2 am \/2Im an \/2Pm a a . 
V2/ + 1 I \fdi \fd, \fE \fd. 



it holds that, for I large, it holds approximately that Ui ~ U(S 21 ), i.e. U\ it is asymptotically 
distributed on the unit sphere of M 2 ' +1 . Under these conditions, it is simple to show that HFE =^ 
HFG, i.e. 



( lim E [C l - l) 2 = ol { Tl{x) ^ N(0, 1) , as I 

Indeed, it is sufficient to note that, as before 

Ti{x) _ Ti imp y/kitaiQ 
y/Var{Ti) ~ ^(21 + 1)Q ~ VQ ' 

which we can write as 

am 



DC 




Now, as I — ► oo 



because the left hand side can be viewed as the marginal distribution for a uniform law on a 
sphere of growing dimension; the latter is asymptotically Gaussian, as a consequence of Poincare 
Lemma (see [10j). We do not investigate this issue more fully here, and we leave for future 
research the determination of general conditions such that (compare with (|29p ) 

the law of Ui and U(S 21 ) are asymptotically close as I — > oo . (30) 
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Obviously, for all fields such that ([30]) holds (i.e. those that are asymptotically completely 
random, to mimic the terminology of [H [9]), by the same argument as before we have that 

- 1 HtpW^ (0 ' 1) }' 

To conclude this work, we wish to provide an example were the HFE and HFG property are 
indeed not equivalent. Consider the (anisotropic) field 




K x ) = ^2iimXim{x) , where & m = I 

7™ V 



lm 

and the random variables verifies the assumption 



£l , for m = 
, otherwise 



= 0, ^Etf < oo and E$f < 



x 



Note that, in the definition of h (x), the sum is not taken with respect to I. The field can be 
made isotropic by taking a random rotation T{x) = h(gx), where g is a random, uniformly 
distributed element of SO(3). We have as usual T(x) = Y^i Yl m =-i a lmXlm( x ) , where 



law / \/- l aw I ^ v I \ r 

aim = 2^ D m'm{9)t,lm' = V 2T+T ' 
rra'=— I 

and where {D l (g)} denotes the well-known Wigner representation matrices for SO(3), and the 
first identity in law is discussed for instance in [2J, Note that 

l . i 

4tt 



21 + - 

m=—l m=—l 

as expected, because the sample angular power spectrum is invariant to rotations. Of course in 
this case we do not have ergodicity in general, i.e. it may happen that 



u |2 c 
l^m=-l \ a lm\ t 



EE l m =-i\ai m \ 2 Eti 
and indeed for general sequences {£/;} 



2 



2 



However, in the special case where 

e~ l with probability | 
— e - ' with probability i 



6 

, • 1 2 



we obtain easily that E j — 1 j =0, while asymptotic Gaussianity fails. Hence, we have con- 
structed an example where the HFE property holds but the HFG property does not. Note that 
the support of the vector is concentrated on a small subset of the sphere S 21 ; heuristically, 
this is what prevents Poincare-like arguments to go through. 
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